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We study the mutual interactions of simple, parallel polymers and nanotubes, and develop a 
scheme to include the van der Waals interactions in the framework of density functional theory 
(DFT) for these molecules at intermediate to long-range separations. We primarily focus on the 
polymers polyethylene, isotactic polypropylene, and isotactic polyvinylchloride, but our approach 
applies more generally to all simple polymers and nanotubes. From first-principle DFT calculations 
we extract the electron density of the polymers and their static electric response. We derive explicit 
expressions for the van der Waals interaction energy under simple symmetry assumptions. 
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I. INTRODUCTION 

Developing an ab initio understanding for soft-matter problems, such as found both in the 
biosciences and in carbon-based nanotechnology, represents an important challenge. Density func- 
tional theory (DFT) provides, in principle, the necessary materials-theory account but traditional 
implementations of DFT are insufficient for soft matter 0, Q ■ Traditional DFT implementations 
provide an accurate account of the intra-molecular (intra-material) binding in regions with a high 
electron density. However, the typical soft-matter problem, like the binding of graphitic and lay- 
ered materials |3|, of polymer matrices |5|, and of polyaromatic dimers |a, is defined by 
regions of sparse electron density. Here traditional DFT implementations provide no consistent 
account [E IS El EOl because they lack a description of the van der Waals (vdW) or dispersion 
forces. This shortfall of traditional DFT has motivated a series of development works by us 
and by others [H O El 13 that seeks a consistent extension of traditional DFT (accurate for 
intra-molecular binding) with electron-density based calculations of the van der Waals interactions 
(for the intermolecular binding) . We thereby extend traditional DFT in a consistent approach that 
avoids double counting of the electron interaction effects. 

The exchange-correlation energy E xc [n] is a functional of the electron density n and represents 
the central concept in all DFT approximations. It represents the total electron-electron interaction 
energy which is not included in the mean-field approximation described by the electron Poisson po- 
tential 3>(r). Together with the external potential V ex t, the functional -E^cM completely determines 
the general DFT calculational scheme for the total energy 

E = T [n}+ J dvn{v) [$(r)/2 + V ext ] + E xc [n] (1) 

of an interacting soft/sparse or traditional- materials system. The term To[n] in is the kinetic 
energy of a non-interacting system at density n and is in DFT straight forward to calculate as 
a single-particle Schrodinger equation problem specific to the materials system, as discussed, for 
example, in Refs. [f5lll6|]. The functional E xc [n], however, is universal and it is this functional that 
we seek to improve beyond the traditional local-density and generalized gradient approximations in 
the investigations for a new van- der- Waals-density functional (vdW-DF) |2ll3ll6lll7|. 
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In this paper we illustrate the nature of such a new vdW-DF by ab initio and analytical calcula- 
tions of the intermediate-range interactions between parallel (unbranched) polymers and nanotubes. 
Specifically, we investigate the interactions between polyethylene (PE), isotactic polypropylene (PP), 
isotactic polyvinylchloridc (PVC), and single- walled carbon nanotubes (CNT). We first determine 
the structural geometry of the isolated molecules by self-consistent DFT calculations using the gen- 
eralized gradient approximation (GGA), allowing the atoms to relax to the energetically optimal 
positions. In the optimization the atoms are moved according to the Hellmann-Feynman forces in 
the GGA; details of these calculations will be published separately [Isj. To this state-of-the-art 
description of the individual molecular structure we here add a charge-density based description of 
the intermolecular binding. 

The elongated molecules represent an important class of soft-matter interaction problems be- 
cause they have very anisotropic dielectrical response functions, and as such constitute a special 
problem. This interaction is naturally separated in three regimes: (i) at binding distances where the 
van der Waals attraction and kinetic repulsion compete, (ii) the intermediate regime where the local 
electron density variation is important but where the electron overlap can be neglected, and (iii) 
the asymptotic regime where the van der Waals interaction is given by the overall (as opposed to 
the locally varying) susceptibility of the molecules and where traditional van der Waals calculations 
[T^. |20| become feasible. We focus on describing the transition of the molecular van der Waals 
interactions from intermediate separations and out to asymptotic separations. 

The emphasis on intermediate molecule separations d permits us to use a simpler approximative 
treatment of the anisotropic dielectric response, which we determine upon averaging the electron 
density along the polymer and nanotube axis. Our aim is to find general trends in polymer-polymer 
and nanotube-nanotube interactions, and rather than aiming at numerical, very accurate results 
we here develop analytical, approximate expressions for use in further analyzing the system. More 
accurate numerical studies, with improved calculations of the molecular dielectrical response, are 
deferred to a forthcoming paper. 



II. THE INTERMOLECULAR VDW BINDING 



The vdW interactions involve fluctuating (dipole and multipole) interactions which are nat- 
urally expressed in terms of the bare (local-field) and external-field susceptibility tensors, xo and 
Xcff- Starting from the adiabatic connection formula [2il . l22T.l23j the dipole-dipole approximation of 
the exchange-correlation energy of two well separated objects 1 and 2 can be written as |24| 



E%™[n] = E vdw [n] ~ - jH J d 3 n J d 3 r[ J d 3 r 2 J d\' 2 



^ A( r i' r i;«) T i3( r l. r a)x^, 3 (r a y 2 ;«)I^(r2,r 1 ). (2) 



i.b.c.d 



Here we allow for anisotropic effective susceptibility tensors of the objects. The anisotropy must 
necessarily be taken into account for such extended objects as the straight polymers and nanotubes 
discussed here. The interaction tensor is given by = — \7iVj\rj — r.j | — 1 . A standard frequency 
integration has been turned into an integral on the imaginary axis, with u = —iuu and u> the physical 
frequency. Hartree atomic units are used throughout, unless otherwise noted. 

In this study we seek only to explore the formal nature of the interactions between parallel 
polymers and nanotubes. We therefore make use of a simplified but electron-density based model de- 
scription of the bare molecular susceptibility xo which is assumed diagonal. This model assumption 
implies that the anisotropic dielectric response Xefl (that characterizes the external-field response 
of polymers/nanotubes) is taken to arrive exclusively by the local-field effects that are specified by 
this diagonal xo- 

For the diagonal elements of xo we use the approach of Refs. 0,|2^|, that is, a slightly modified 
version of the effective susceptibility used in Ref. 
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Xo (t,t>;u)=6(t-t')^L. (3) 

IT + Uq 

The valence electron density n(r) is obtained directly from the above-mentioned self-consistent DFT 
calculations. The cut-off frequency uq was introduced in Eq. © to avoid an unphysical divergence 
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FIG. 1: Valence charge density n(s) of polyethylene (PE), polyvinylchloride (PVC), and the _R( g0 )-nanotube 
(CNT), averaged along the length of the molecules. Contour line spacing is 0.15 eA -3 and only densities equal 
to or above this value are shown. The valence charge density plot of polypropylene (PP) looks qualitatively 
similar to that of PVC. 



in the static limit (u = 0). The cut-off is determined from a separate set of self-consistent DFT 
calculations that provide results for the static polarization by imposing a finite external electrical 
field pl.l2f)| perpendicular to the polymer and nanotube axis. 

Figure n shows the electron density (averaged along the polymer and nanotube axis) as calcu- 
lated by traditional DFT for the molecules PE, PVC, and CNT. The averaging along the length is 
done in the same spirit as our calculations for the vdW interactions in layered materials UEM 
The averaging provides a qualitative picture of the electron-density variation that effectively deter- 
mines the intermolecular vdW interaction at the relevant intermediate molecular separations. For 
simplicity we choose to base the set of numerical and analytical evaluations of the vdW correction 
© on dielectric-response functions specified by these averaged electron densities. 

Figure n a l so suggests that this effective electron distribution has an approximate cylindrical 
rotational symmetry which we utilize for further simplifications in our present model study. We term 
such molecules 'cylindrical' and these remain the focus of this paper because they offer in several 
cases comparison against analytical expansions and evaluations. We stress that for the question 
of the mutual vdW interactions the most relevant measure of cylindrical symmetry is exactly the 
angle-dependence of the static external field response. In the set of investigated molecules this 
angle-dependence is small, as shown for PE in the inset of Figure EJa) from direct evaluation in 
traditional DFT. 

An element of charge at position s (specified in two dimensions relative to the molecule center- 
of-mass line) on the molecule responds to the local electric field according to Xo( s ; u) given by J3J. 
The local field is the applied electric field E app ii e( j(u) (which arises from charge fluctuations at the 
other molecule) screened by the presence of other elements of charge in the neighborhood. The local 
electric field is described via the local electric potential cf>(s; u) and given by charge conservation 

-V-{(l + 47rxo)V0} = O . (4) 

The effective susceptibility tensor Xcffj which appears in the expression for the interaction energy 
@ , is related to the local susceptibility and the local electric potential by 

Xes (s; w)E applicd (u) = -Xo(s; u)V0(s; u) . (5) 

Due to the assumed translational invariance along the length of the molecule we can in general 
solve (QJ numerically as a two-dimensional differential equation in </>(s; u) with the frequency u as 
a parameter, and thus find Xefi( s ]u) from (jSJ. The spatial integral a{u) = L J d 2 sxcfi (s; u) taken 
at u — gives us the molecular, macroscopic susceptibility tensor as a function of uq. Comparing 
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FIG. 2: The van der Waals part of the interaction energy for two parallel polymers at separation d. Panel 
(a) shows E v dw/L numerically evaluated for sets of length-averaged polymers keeping their full angular 
structure. The squares in the inset show the DFT calculated susceptibility of PE for different angles. The 
full line is a fit based on the rotational properties of the susceptibility tensor. Panel (b) shows the results for 
the vdW interaction energy of PVC in different approximations. The thick lines correspond to evaluations 
keeping the full angular structure for the length-averaged electron density, and the thin lines correspond to 
evaluations of the radially averaged density. Shown are expansions to 5th order in d" 1 (dotted lines), to 7th 
order (dashed lines), and the full calculation without expansion (solid line). In both panel (a) and (b) the 
PVC molecules are oriented as shown for separation d = 10 A in the inset of (b). The PP molecules are 
oriented similarly to PVC, and the PE molecules are oriented with the long axis towards each other. 



a(u = 0) to the static DFT-calculated molecular susceptibility we determine the frequency constant 

In polymers our approach [J, |25| yields a value of the static polarization which varies with the 
angle at which the electric field is applied. However, for the simple polymers the variation is small, 
as shown for PE in the inset of Figure |2f a) . For PP and PVC the susceptibility is almost constant 
with respect to rotation, with a pp = 0.97 ±0.01 eA 2 /V and a PVC = 1.01 ±0.01 eA 2 /V. The values 
uq found for the angular means of the static polarizability are 0.429 Ha for PE, 0.507 Ha for PP, 
and 0.482 Ha for PVC. In nanotubes uq depends on the nanotube radius. However, the variation 
is small and the macroscopic susceptibility is not very sensitive to this variation. For CNT we thus 
use the large-nanotube value |2f| for all nanotubes, uq = 0.30 Ha. 

In the expression for the van der Waals energy the two spatial integrals along the (density 
averaged) molecules can be performed analytically [4|, but leave a cumbersome integral over the 
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remaining spatial variables and the frequency. These integrals can be numerically solved, and the 
results for pairs of parallel polymers are shown in Figure |5Ja) for a specific relative orientation of 
the polymers. The results for other relative orientations will be reported in Ref. [l£l |. 



III. INTERACTION ENERGY IN CYLINDRICAL SYMMETRY 



Whereas it is possible to numerically evaluate E vt xw for a general pair of parallel molecules, 
such results provide only limited insight into the vdW interaction mechanism. To interpret such 
numerical evaluations we provide instead analytical results that are possible to derive with our choice 
of simplified local susceptibility function for cylindrical molecules. When cylindrical symmetry 
is imposed, the local electric potential tj> factorizes into <fi(s;u) — —\E app u e a(u)\W(s;u) cos(9, where 
9 measures the angle from the direction of the applied field. Solving 0} then reduces to the one- 
dimensional problem of finding W(s; u) with u as a parameter. In Ref. 18] the more complex 
problem of non-cylindrical polymers will be treated in greater detail. 

For cylindrical molecules XeS is diagonal in the cylindrical representation (s, 8, z) |25j, and the 
interaction energy reduces to 

/>oo 

dsisi j ds 2 s 2 ^2 xf f f(si; w )Xcff( s 2;u)G / 3 7 (s 1 ,s 2 ) (6) 
where the geometry factors are 

/•2tt />2tt />oo 

G 01 ( Sl ,s 2 ) = de x \ d9 2 d(z 2 -z 1 )[T^(s 1 ,e 1 ,z 1 ,s 2 ,6 2 ,z 2 )] . (7) 

JO JO J-oo 

The integral over z 2 — z\ [4| and one of the angles in G^ 7 (si, s 2 ) can be carried out analytically, but 
still leads to intractable expressions in the remaining angle, si, s 2l and d, and the full presentation 
is deferred. 

Our further treatment of the interaction energy expression differs depending on the nature 
of the molecules considered. On the one hand the polymers considered here are thin (radius ~ 
2 A) compared to the typical nanotube (radius w 10 A) and compared to the typical binding 
separation (« 3-4 A plus radii of molecules). This means that a large-separation approximation, 
based on the center-of-mass to ccntcr-of-mass distance d, is a better approximation for the polymers 
than for nanotubes. On the other hand, nanotubes are about the most cylindrically symmetric 
extended molecules that exist, and the major part of the electron density is concentrated in a thin 
cylindrical shell, as shown in Figure Q It is therefore reasonable to approximate the nanotube 
electron charge density by a weighted radial delta-function. We further approximate the nanotube 
effective susceptibility tensor by assuming xts = xts- Below, we pursue these approximations for 
polymers and for nanotubes. We also show that the two sets of derivations are consistent under 
relevant conditions. 



L Jo 2tt L 



IV. POLYMER INTERACTION 



As argued above, we can explore the polymer interactions in an expansion in powers of d 1 . 
Expanding the geometry factors, writing — 2J 4 G '^'d , the first three even-index terms are 

G (o) = ^ N (8) 
128 v ' 

G(2)(S1 ' S2) = 2 -^((N + 2M)sf + (N + 2Mfsl) (9) 
G(4)(S1 ' S2) = 3 -^((N + m4 + 3L.sl S l + (N + AM) T 4) (10) 

with the tensors 

/ 19 19 10 \ / 12 12 8 \ / 99 73 52 \ 

TV = 19 19 10 , M = 7 7 2, L = 73 59 28 . (11) 
\ 10 10 12 J \ 5 5 6 / \ 52 28 48 / 
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For cylindrical molecules all terms odd in i for symmetry reasons vanish after integration in 
© , the odd-index terms being a measure of the molecular charge anisotropy. The polymers treated 
here in reality are cylindrical symmetric to a varying degree. We find that for PE interactions the 
numerical values of the first few odd- index contributions to the energy are indeed small 0] . 

To lowest order in dr x the interaction energy becomes 



8,-y 



9tt 2 



256 d 5 ^ 



Nf} y J 07 (12) 



where we have defined the frequency integral 



J/3 7 = / d«Xcff,l( U )Xcff. 2 ( W ) ( 13 ) 





given by the spatial integral of xfg 

= 2^1 / d3 ^cff(«; «) = J o °° daa&fa u) . (14) 

In the approximation \ls(u) — xts ( u )> use d below for nanotubes, we find from l|12(l the asymptotic 
contribution to the interaction energy 



E {0) 9ir 2 



L 256d 5 



(4 ■ 19 J ss + 2 ■ 10J SZ + 2 ■ 10J ZS + 12J ZZ ) (15) 



Using @ and (|10l) and the corresponding generalization of (|13fl and 114fl the next order corrections 
can be similarly expressed. 

In Figure 0;b) the cT 1 -expansion is illustrated along with a numerical calculation that keeps 
the full angular dependence not assuming any rotational invariance (solid line) . The energy is for the 
PVC-PVC interaction in the relative arrangement shown in the inset. This is the orientation affected 
the most by the anisotropy of the polymer. The thin dotted line is the lowest order contribution 
(|12l) of the expansion. As the next non- vanishing order, d~ 7 , is added the description is improved 
(thin dashed line). We checked that adding the 9th order did improve the description of the energy 
slightly, but not significantly at the separations considered here. As a comparison, we also show the 
d^ 1 -expansion of the full angular calculation: the thick dotted line shows the d~ 5 term, and the 
thick dashed line includes all orders up to dr 1 . We note that using the assumption of cylindrical 
symmetry of PVC does worsen the quality of the evaluation somewhat, but already including a few 
terms in the <i _1 -expansion improves the approximations. 



V. NANOTUBE INTERACTION ENERGY 

For nanotubes a useful simplification of © and J7J) can be obtained by a different approxima- 
tion. Rather than using the d - ^expansion, we stress that the nanotube electron density is relatively 
well described as a weighted radial delta- function located on the nanotube radius R [25j . This leads 

to the approximation xf ff (s; u) « Xeg( u )H s ~ which decouples the susceptibility factors from 

the spatial integration, and we find 



-EvdW 

L 

f3,~f=s,6,z 



E / f=Ai(«)3^(«)G/»r(i*i,ik) 

l=s,8,z J ° 

= 4 E JfhGM{RuRa). (16) 
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This is the expression for the interaction of two parallel nanotubes of radii R\ and R 2 . The integrals 
in Gfj 1 {R\, R 2 ) yield, in general, tedious expressions with one remaining angular integral. However, 
with the introduction of another, numerically small approximation of the effective susceptibility 
tensor we can carry the analysis of E v( iw further with simple expressions. 

Assuming that the averaged effective radial and tangential susceptibilities on each molecule 
are identical, xtg ( u ) = xlg( u )> tnen Jss = Js6 = Je s = Jee, Jsz = Jez, and J zs = J z8 . This is a less 
restrictive approximation than the already modest approximation Xeff( s > u ) = Xcff( s ! u ) imposed in 
Refs. The van der Waals energy per nanotube length then becomes Q 



' ' ' G ss + G s g + Gg s + Gee ) + Jsz ( G sz + G 



+ Jzs(G zs + G z o)+J zz G zz } (17) 

where Gp 7 is implicitly taken in nanotube radii R\ and R 2 . 

In Eq. (fT7|) all the sums of G-terms are of the same form, G ss + G s g + Gg s + Gee = 19G to t, 
G sz + G Sz = G zs + G zS = 5G to t and G zz = 3G to t where 



G tot (R u R 2 ) = ^- [ dOx [ d9 2 ( (R 2 cos 6 2 + d - R x cos 9 1 ) 2 + {R 2 sin 9 2 -R x sind 
^ Jo Jo 
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r w/2 

/ 2F1 (5/2, 5/2; 1; {{R 2 - i?x) 2 + 4R X R 2 sin 2 £}/d 2 ) . (18) 

Jo 



The second line is a significant simplification which expresses the geometry variation in terms of the 
hypergeometric function p F q integrated over the variable £ = (0 2 — #i)/2. 

For i?i = R 2 = R also the remaining integral can be carried out, again in terms of a 
hypergeometric function |2y|. We can thus factorize E v dw into one frequency integral Jtot = 
19 J ss + 5 J zs + 5J SZ + 3J ZZ and one spatial integral G to t(R, R) 

— z — = — jr—Jtot Gtot{R, R) 

= 3^ 2 (1/2, 5/2, 5/2; 1, 1; 4i? 2 /d 2 ) . (19) 

The only numerical integration needed is the frequency integral in J tot . This simplification previously 
allowed us to extract the graphite-graphite limit of two large-radius nanotubes at intermediate 
separation 0. It also allows us to interpret the above-mentioned d _1 -expansion of cylindrical 
molecules: By expanding (|19fl to the lowest non- vanishing order in d~ l we recover ()15fl . 



VI. CONCLUSIONS 



We have here presented a general approach based on first-principles electron density calcu- 
lations for computing the intermediate to long-range interactions of parallel, geometrically simple 
polymers and nanotubes. We have applied the scheme to polyethylene, isotactic polypropylene, and 
isotactic polyvinylchloride, and the interaction energy for these molecules has been evaluated. We 
argue that the polymer electron charge density can be approximately described as rotationally in- 
variant, which allows us to derive a number of simple explicit expressions for qualitative insight. For 
parallel nanotubes we derive expressions for the interactions in terms of hypergeometric functions. 
Further, the asymptotic expansion for the nanotubes is seen to be consistent with the expression 
obtained for the cylindrical polymers. 
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